A model for QCD ground state with magnetic disorder 
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We explore an ansatz for the QCD vacuum in the Coulomb gauge that describes gauge field 
fluctuations in presence of a weakly-interacting gas of abelian monopoles. Such magnetic disorder 
leads to long-range correlations which are manifested through the area law for the Wilson loop. 
In particular we focus on the role of the residual monopole-monopole interactions in providing 
the mechanism for suppression of the gluon propagator at low-momenta which also leads to low 
momentum enhancement in the ghost propagator. 
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I. INTRODUCTION 

A condensate of magnetic monopoles screens chromo- 
electric field between fundamental color charges and leads 
to formation of flux tubes between static quarks. Con- 
sequently condensation of magnetic degrees of freedom 
in the QCD vacuum has been proposed as the mech- 
anism underlying color confinement [1-4] and recently 
strong evidence for magnetic dominance has emerged 
from lattice gauge simulations [5-8]. In QCD monopoles 
emerge as solutions of classical equations of motion [10], 
albeit having infinite energy, the underlying ultraviolet 
singularity is expected to be regularized by quantum 
fluctuations [11, 12] and in any case should not play a 
role in the low energy domain. Monopoles create vor- 
tices [5, 13, 14] and all together lead to percolation of 
monopole-antimonopole chains [15]. Even though there 
is ample evidence for presence of magnetic domains it is 
still an open issue which of the many possible monopolc- 
vortex geometries dominates the QCD vacuum [8]. It 
should be noted however, that vortex-like configurations 
of long monopole chains arc needed to achieve the confin- 
ing scenario based between charges associated with the 
center of the gauge group as expected for QCD [9]. 

Recently lattice simulations in both Landau [16-24] 
and Coulomb gauge [25-27] have significantly advanced 
our knowledge on the infrared (IR) properties of QCD 
Green's functions. And together with studies in the con- 
tinuum [28-32, 34-38, 51] have reinvigorated discussion 
on the role of Green's functions in probing the long-range 
properties of the QCD vacuum [39, 40]. For example, in 
the Gribov-Zwanziger (GZ) scenario [41-44], confinement 
is related to presence of large field configurations near the 
boundary of the Gribov region, the Gribov horizon. The 
Gribov region is defined as the domain in the gauge field 
space that satisfies a particular gauge condition. Within 
the Gribov region the Faddeev-Popov (FP) operator is 
positive and it vanishes on the horizon. Thus in the 
GZ scenario one expects that the (vev) of the inverse 



of the FP operator, referred to as the ghost propagator, 
is IR enhanced. On the other hand, the gluon propa- 
gator, which represents propagation of color charges is 
expected to be suppressed. The IR enhancement of the 
ghost propagator and suppression of the gluon propaga- 
tor have unambiguously been seen in lattice simulations 
in both Landau and Coulomb gauges. This is to be con- 
trasted with some of the phenomenological models were 
it is assumed that the quark-antiquark potential origi- 
nates from the nonperturbative gluon exchange. Thus a 
confining quark-antiquark potential would necessitate an 
IR enhanced gluon propagator [45, 46]. 

In this paper we explore the possibility that the low- 
momentum suppression of the gluon propagator may not 
necessarily be related to confinement but to color screen- 
ing. In the latter case IR suppression of the gluon prop- 
agator implies that colored, physical gluons do not prop- 
agate and as follows from the GZ conjecture it would 
not be the role of the gluon propagator but of the IR 
enhanced ghost propagator to carry the distinct signa- 
tures of confinement. It is a known phenomena that a 
dynamically generated gluon mass and IR suppression of 
the gluon propagator can emerge as a result of screening 
e.g induced by condensation of magnetic domains, (vor- 
tices, monopoles) [47, 48]. It can be well illustrated in 
models [8], i.e the Fradkin-Shcnker model [49] where the 
confined and Higgs phase are smoothly connected and 
there is also a smooth transition in the gauge propaga- 
tors [50]. Similarly in the Landau gauge, solutions of 
QCD Dyson-Schwinger equations that display screening 
behavior have recently been studied in [51]. 

In the following section we examine the effect of con- 
densation of magnetic monopoles in the Coulomb gauge. 
In particular we study an ansatz for the vacuum wave 
functional that contains a weakly interacting gas of 
monopoles. The recently studies of various models for 
the vacuum wave function [34-38, 52-55] were shown to 
be quite successful phcnomenologically. In particular an- 
alytical calculations can be done with an ansatz that 
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describes the vacuum in terms of gaussian fluctuations 
of the transverse, vector potential around the zero-field 
configurations. Since such a simple vacuum does not ex- 
plicitly contain magnetic degrees of freedom it leads to 
a perimeter law falloff of the vev of a Wilson loop. Nev- 
ertheless onset of confinement could be seen through the 
1R enhancement of the ghost propagator and the non- 
abelian Coulomb potential between static charges. This 
is because gaussian ansatz can support large fluctuations 
of the fields, which may be nearing the boundary of the 
Gribov region. We begin with the construction of our 
vacuum ansatz followed by a discussion of the gluon and 
ghost propagators as well as the Wilson loop. A sum- 
mary and outlook are given in Section III. 



II. MODEL FOR THE MONOPOLE 
DOMINATED VACUUM 

In the following we consider SU(2) Coulomb gauge 
QCD, in the Shrodinger picture represented by the trans- 
verse field variables, V • A a = 0, where a = 1 • • • Nq — 1 = 
3 is the color index. The gaussian approximation to the 
vacuum discussed earlier is given by a wave function in 
the from 



an abelian monopole centered at c is given by [57] ( x n = 

TX * X j X \^ - — X TX X 



a la (x, a) = qa l (x - c, n)w a , 



where 



a l {x — c, n) = 



[{x — c) x n]' 1 



4tt \x - c| [\x - c| - (x n - c„)] ' 



(4) 
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a = (q, c, n) denotes collectively the monopole coordi- 
nates, q ± 1 is the monopole charge in units of the mag- 
netic charge <?, c represents its location and n is a unit 
vector that defines the orientation of the (straight) Dirac 
string. In the following we will neglect fluctuations in 
the relative orientation of the monopole color orienta- 
tions, i.e in Eq. (1) we set w to be the common orien- 
tation for all monopoles. This restriction can be easily 
removed, however, analytical calculations, even in the 
weak coupling limit that we discuss below, would not 
be possible. The abelian component of the gauge field 
A L is then assumed to fluctuate over a background of 
monopole-antimonopole gas. For N monopoles we thus 
write {a = (a\, ■ ■ ■ ajv)), 

* lIv [^4 l ] a e - k I dxdyA L (x,a)uj L (x-y)A L (y,a) ^ ^ 



where 



V%[A 



t oc e 



■| / dxdyA T (x)u> T (x-y)A T (y) 
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We will use the T(L) subscript to represent components 
in the color algebra transverse (longitudinal) to a chosen 
direction, w = w a , (w 2 = I), i.e. Aip = A ia - (w b A tb )w a , 
A l L = A la w a . The meaning of w and the role of fields 
along w will be discussed below. The 2-point func- 
tion (hereafter referred to as the gluon propagator) of 
color-transverse gluons is then given by (V is the three- 
dimensional volume ) 



A L (x, a) = A L (x) - ^ a ( x > 



(7) 



Finally, after summing over the monopole coordinates, 
averaging over their (common) color orientation and 
summing over the N configurations we obtain the ansatz 
for the vacuum wave functional given by, 
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Z[0] 



dw 



Da^* m (a)|^;«[A L ]^[A T ]| 2 



N=0 



, (8) 

V (A T (k)A T (-k)) = S T (k)5 T (w)D T (k), (2) Thc intcgration mcasurc over monopole coordinates is 



where 6t(ti) = 5jl(ri) = 6ij — n % n^ /n 2 is the transverse 
projector in three dimensions and Dr(k) = l/(2wr(fc)) 
where wr(fc) is the Fourier transform of ujt(x) that ap- 
pears in Eq. (I). Our ansatz wave function absorbs the 
Coulomb gauge functional measure [56] i.e. 

J DAJ[A]^ QCD \A) 2 -+ J DA{* Ansatz \A) 2 = 1. (3) 

With the gaussian wave function to (approximately) in 
order to restrict field configurations to thc inside of the 
Gribov volume where the FP operator J is positive, 
ujr(k) should be taken such that it diverges in the IR 
limit i.e. u T (k -> 0) cx k 1+a °- T with a G , T < -1. For 
large momenta, the choice ur(fc) — > k makes the vacuum 
wave function match that of the free theory. 

The separation of transverse and longitudinal, (with 
respect to w) components is motivated by the assumed 
dominance of abelian monopoles. Thc classical field of 



given by Da = II 



- tin 

i=l 



The 



dc,(dn 4 /4^)(I/2E, ?l=± i 
distribution of monopoles is specified by the wave func- 
tion $ m ({a}). We first discuss the non- interacting ap- 
proximation, 



*m(«) 



const., 



(9) 



with p being the density of monopole pairs. The gluon 
propagator in this case is simply given by, 



V-V'W^ (-*)) = S^(k)5 ab D(k), 
D(k)= 2 -D T (k)+ l -D L {k) 



(10) 



with 



D L {k) 



1 



2uj L {k) 
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2u L (k) 
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dn 



a l (k,n)a l (-k,n)] 



dx- 



(11) 
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The monopole contribution is singular in the limit e — > 
due to the collinear singularity associated with momen- 
tum component along the Dirac string [57, 58]. Even 
if this divergence was to be regularized, for example 
by combining monopole-antimonopole paris into closed 
chains, the contribution remains strongly enhanced in 
the IR due to the 1/k 4 behavior which originates from 
the long-range, Coulomb, monopole field. Thus a sim- 
ple model with non-interacting monopoles (same result 
is obtained in the case of vortices) cannot be adequate 
since it gives as strongly enhanced gluon propagator that 
is inconsistent with all lattice results. 

The IR suppression of gluon propagator must therefore 
originate from screening by the interacting monopoles. 
To this extent we introduce an effective interaction which 
is repulsive (attractive) between monopole and (anti) 
monopole, respectively, 

Vij = V(ai,aj) = qiqjS 2 (ni -rij)V{ci -Cj,rii), (12) 

and replace by the corresponding partition function 

JV 

* m (a) = ^e" (13) 
We choose the potential V in Eq. (12) in the form 

V{k,n) = J dxV(x,n)e^ = ^p^- (14) 

As will be seen below, the 1/k 2 , (kz = n ■ k) term will 
be responsible for screening in the direction along the 
string white the critical exponent 7 will control the IR 
behavior of the propagator, and M will related to the 
inverse Debye length. 

In calculation of matrix elements, (O) = 
J DAO[A}\Hf[A]\ 2 with the wave function given by 
Eq. (8) the summation over the number of monopoles is 
done with the help of an auxiliary sine-Gordon field [3], 
4>{x,n). In particular for a generating functional 

Z[J] = J DAef dx ^ rL ^ A L^ +J T^ A T( x )]\-q/[A]\ 2 , (15) 

one finds 

Z[J]=Z T [J T ] j D0e _,Sl[ * , " 7l] , (16) 

with Zt[Jt] — exp(| J dxdy,JT(x)uj^, 1 (x — y)JT{y)) and 

= / dxdn [4>{x,n)d 2 L {d 2 Y + ^{x,n) 

+A D (4+27) [1 - cos^z, n) - v L (x, n))}] , (17) 



We defined the Debye screening length \d = 
(M 1+27 p)~ 1// ( 4+27 '. In terms of the generating func- 
tional the gluon propagator is given by d 2 hiZ[J = 
®]/dj(k)dj(-k)- I n th° mean-field approximation, which 
is valid in the limit Xd > p~ 1//3 , where p -1 / 3 is the av- 
erage separation between the monopoles, the functional 
integral over </> can be computed analytically using the 
saddle point approximation. At this level it is clear that 
interaction between monopole pairs screens their charges 
and thus regularizes the IR divergent part of the gluon 
propagator, leading to a propagator which is of the form 
given by Eq. (11) but with the last term on the right 
hand side replaced by 



where 



v L (x,n) 



9 2 P 



J dx{\ - 



2) fc 27 

X ' k 2 n k 2 + 2 ^ + A" 4 " 27 

(1 + k%) arctan(A;f>) - 
P 



k 2 

D 



,(19) 



where ko = kXp. In [27] lattice simulation of the 
Coulomb propagator in D — 4 was performed and shown 
to be well approximated by the Gribov formula, D(k) = 
k/(2 v / k i + M 4 ) with M = 0.890 GeV. In [53-55] it was 
shown that the D = 3 YM action in Landau gauge may 
be a good approximation to D = 4 Coulomb gauge wave 
function, and recently the corresponding Coulomb prop- 
agator was evaluated and shown to compare favorably 
with the exact D = 4 propagator [52]. For compari- 
son we then use results obtained with the D = 3 YM 
action since it allows to separate the non-magnetic and 
magnetic contributions. In Fig. 1 circles represent the 
Coulomb propagator from [52] with vortices removed. 
We use this propagator to fix and lot under a simpli- 
fying assumption Wi,(fc) = u>T(k) which we parametrize 
as (k/m) a /2({k/m) 2 + 6 2 ) 1 /2+a/2 so t nat a t low momen- 
tum ujl.t ~ k a and in the UV, ujl,t ~ The fit 
of lu to the propagator from [52] with vortices removed 
yields, a = 1.21, b = 0.29, m = 3.88 GeV. The re- 
sult of the fit is shown by the dashed line. We then fix 
the three parameters that describe our propagator (c.f. 
Eqs. (11), (19)), L = g 2 p\%, \ D and 7 by fitting the full 
gluon propagator from [52] (squares in Fig. 1) and ob- 
tain X D = 1.51 GeV" 1 , L = 6.46 GeV" 1 and 7 = 0.32. 
The result is shown by the solid line. And the condition 
of applicability of the mean-field approximation requires 
week monopole-monopole interaction ( i.e. strong chro- 
moelectric coupling, e = kx/g) with 




< 2. 



dyJ z L (y)a l (y - x,n). 



(18) 



(20) 



As discussed in Sec. I, and shown above, while the IR 
suppression of the gluon might be the result of screening 
of magnetic charges, IR enhancement of the ghost prop- 
agator comes from YM field distribution near the Gribov 
horizon. We postpone a detailed numerical study of the 
ghost propagator and just notice that because monopoles 
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FIG. 1: Comparison of our gluon propagator with that ob- 
tained from lattice computations [52] 



introduce orientation in color, the diagonal (-Dl) and off- 
diagonal (Z?t) gluon propagators are different, and the 
mean-field relation between the ghost and gluon propa- 
gators becomes more complicated. In particular, defining 
the longitudinal and transverse ghost form factors di(fe) 
and dx{k) by 



d L (k) _ 
k 2 ~ 

d T (k) 
k 2 



1 J DAw a w h M- 1 [A]\^ w [A]\ 2 
v -i J DA S rM M - 1 [A}\^ w [A}\ 2 , 



(21) 



where M~ l [A] = Af- 1 [A](fc, a, -k, b) = -e(V • £>[A]) _1 
is the inverse of the Faddeev-Poppov operator. In the 
mean-field approximation one obtains, 

dl 1 (k)=e- 1 -I[D T ,d T } 



d T \k) 
where [59] 



I[D L ,d T }--I[D T ,d L ], (22) 



nD,d] = fJd P D( P )^± d ( P - k) . (23) 

In the case D L = D T = D these give d L = d T = d with 
d given by 



d^ik) =e" 1 -I[D, d] 



(24) 



that has been studied in the past [34-38]. In general, 
after averaging over monopole color directions, 



d(k) = ^d T (k) + l -d L {k). 



(25) 



The IR analyst of Eq. (24) leads to a relation between 
gluon and ghost critical exponents. That is assuming the 



IR behavior, of the form D{k) oc k 21 and d(k) oc k 2S one 
finds, S = —1/4 — 7/2. Thus an IR enhanced ghost prop- 
agator (5 < 0) necessitates a screened and IR suppressed 
(7 > 0) gluon propagator. In the case of Eq. (22), with 
D L{T ){k) oc fc^m an d d L{T ){k) cx k 25L <- T > , respectively, 
one finds St = —1/4—7^/2 and 5l = — 1/4— (7t~ 7l/2). 
In addition one should consider the Coulomb form factor 
and the gap equation. It was found in [38] that the cou- 
pled set of Dyson equations for these functions admitted 
only IR finite solutions. With the addition of monopoles, 
however, a preliminary analysis of Eq. (21) indicates that 
IR critical solutions are possible. 

Finally we comment on the role of monopoles in sup- 
pressing large Wilson loop, 

Wj[C] = / ?Y(vE|Pcxp(ie I dx l A t - a T a )\^) , (26) 
2J + 1 J c 

where T a are the SU(2) color generators in the J-th rep- 
resentation. The integration over A is computed by shift- 
ing Al according to Eq. (7). In the limit of large loops, 
the contribution from the non-monopole component (and 
At) is determined by the gluon propagator, 



In IT, 



[C] ~ - / 
Jc 



C— >oo 



dx'dy^Sfji 



(k)e 



ik-x 



0(R~ 1 D(k ~ i? _1 ), 



(27) 



where R is the perimeter of the loop. An IR suppressed 
gluon propagator with D(k) ~ fc 27 and 7 > leads to 
screening of the Wilson loop, i.e. the loop is dominated 
by shot range correlations and has at most perimeter 
dependence. Thus if the long-range correlations are to 
dominate they must come from the monopole gas and it is 
possible to ignore the contributions from the fluctuating 
field. This leads to, 



W. 



W = 27TT £ 

m=—J 



-S[4>,mn[C]] 



(28) 



where S is given by Eq. (17) and rj[C](c, n) = § c dy l d l (y— 
c,n). Here c is the location of a single monopole. In 
particular for a large loop in x — y plane with perimeter, 
R » \x\, r/ — >• 2-7T sign(c • n). It immediately follows 
that iV-ality zero (J-integer) loops are screened and all 
non-zero Af-ality loops behave equivalently to the loop in 
the fundamental representation, J = 1/2. The Casimir 
scaling presumably comes from the neglected effects of 
fluctuating field. In the weakly-interacting limit (A^ > 
p -1 / 3 ), the path integral in Eq. (28) can also be evaluated 
in the saddle point approximation, this time however the 
saddle point does not correspond to (f> — but centers 
around 77, [3, 60] and is given by the solution of 

(— d ) 2 f— d) 2 ^ 

: ' L - (4+ 2 7 ) (t)(x,n) + sm(<j}(x,n)+mr][C}(x,n)) = 0. 



A 



D 



(29) 

For large loops R >> A in the x — y the equation becomes 
effectively one-dimensional , <j>(x,ri) — > 4>(z n ), z n — z ■ n 
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with 4>{z) interpolating smoothly the discontinuity in r\ 
across the Wlson loop, between —ir and 7r over a distance 
of the order of Xp. Substituting such a saddle point 
solution into the action in Eq. (28) one finds 

i \-4+2 7 

S ^ 2^ ~ R X M^ (30) 

and thus 

lnH/ 1/2 [C] ~ \ DP A[C] (31) 
where A is the area of the loop. 

III. SUMMARY AND OUTLOOK 

Since condensation of magnetic degrees of freedom is 
known to be present in the QCD it becomes essential 
to include magnetic degrees of freedom when construct- 
ing models of the QCD vacuum. In our construction we 
have assumed that these can be represented by a weakly- 
interacting gas of aligned (in color) abelian monopoles. 
Such a state is known to reproduce the confining prop- 
erties for large, fundamental Wilson loops., however it 



does suffer from yielding U{1) N rather then Zjv charge 
dependence. The latter could originate from vortex con- 
figurations of long monopole chains and a construction 
of an ansatz for the corresponding vacuum sate would 
be highly desirable. Here we have shown how magnetic 
monopoles influence the gluon propagator and have ar- 
gued that the IR suppression is the result of screening of 
magnetic charge and not confinement. Suppression of the 
low momentum gluon propagator leads to suppression of 
long-range gluon fluctuations and restricts the field vari- 
ables to the inside of the Gribov region. A self consis- 
tent calculation of gluon and ghost propagators and the 
Coulomb form factor is currently underway. 
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